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By Itiro Tani 


INTRODUCTION - LAMINAR -FLOW AIRFOIIS 


1, In high speed flight conditions, the drag of an airfoil is 
almost exclusively due to skin friction. Therefore, if further reduc- 
tion in drag is desired, it is necessary to delay as much as possible 
the transition from laminar to turbulent flow in the boundary layer along 
the surface, thus decreasing the extent ,of the turbulent boundary layer 
which gives considerable skin friction. As the factors that may affect 
the transition, we will consider the stream t\irbulence, the surface 
rovighness, the surface pressure distribution, and so on. In actual 
flight conditions, however, the effect of txrrbulence seems to be unex- 
pectedly small, so that, so far as smooth surfaces are concerned, there 
remains only the shape of the airfoil section in relation to pressure 
distribution as the most in^jortant factor affecting transition. We call 
a laminar-flow airfoil that airfoil in which the shape of the section 

is suitably designed so as to delay the transition of the boundary 
layer flow. 

2. It is evident that the laminar separation of the boimdary layer 
may cause the transition, as will be mentioned in the appendant part of 
the paper, paragraphs 35-^0* We cannot expect, therefore, to maintain 
laminar flow beyond the separation point. Summarizing the results of 
flight experiments on airfoils hitherto made (refs. 8 to 12), we have 
the conclusion that the observed transition coincides approximately with 
the calculated laminar separation point at small Reynolds nxambers, while 
it moves upstream toward the minimum pressure point as the Reynolds 
number increases. However, no exan 5 )le has ever yet been observed in 
which the transition moves ahead of the minimum pressure point. We 
therefore arrive at the supposition that the laminar-flow airfoil may be 
most simply realized by designing the airfoil In which the minimum pres- 
sure occurs well downstream. 


Kyokaiso no Sen'i o okuraseru Ybkvigata ni tulte." Report of the 
Aeronautical Research Institute, Tokyo In^jerial University, No. 25 O 
(vol. 19, no. 1), Jan. 19**-3- 
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DESIGN OF SYMMETRICAL AIRFOIIS IN WHICH THE MINIMUM 
PRESSURE OCCURS DOWNSTREAM 


3 . Following Professor Moriya (ref. 13) > we write the coordinate 

along the chord in the form x = — (l + cos I) , and assign x = 0, 

I = n to the. leading edge, and x=l, 1=0 to the trailing edge. 

00 

Expressing the ordinate of the mean camber line by M = ^ a^ cos ni, 

0 

and the half -thickness measvired normal to the chord by T = ^ b^ sin nl > 

. ■ 1 

the pressure distribution around the airfoil in the two-dimensional 
potential flow is given by 


£= 1 
<1 
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w 00 I 

q: i sin I + y nsjj (l - cos nl) j ) nbj^ sin nl ? 


+ sin a 4- i a - COS 
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1 + (fAg i eBc)‘ 


iThe ordinates of the upper and lower surfaces are given by M + T 
and M - T, respectively. 
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where a is the angle of attack, p is the pressxire acting on the air- 
foil surface, measured from the static pressure of the undistxirhed stream, 
and q is the dynamic pressvure of the undisturbed streeim. We assume that 

00 

the trailing edge is sharp, so that ) nh^ =0. We limit the range of 

1 

the variables | between 0 and jt, and assign the upper and lower 
parts of the double sign for the upper and lower surfaces, respectively. 
Writing f for the maxlmijm value of M (the maximum csimber) and e for 
the maximum value of 2T (the maximm thickness), we put 


1 - cos nl 




^^s = -2 m nan 


sin nl _ dM 
sin I dx 


^-22- sin n| 

eBs = 2 ^ nbn —— 
1 sin I 


_ ^ ^ cos nl dT 

1 “ sin I dx 


The lift coefficient is given by 


I w 

Cj^ = 2jr ^ sin d - 2 V cos a 


We consider first only the thickness of the airfoil (the cainber of 
the center line will he considered in the next section, paragraphs 8 
to ll) • Namely, we consider the. symmetrical airfoil section set at zero 
angle of attack, with a view to obtaining the minimum pressure well 
downstream. 

4. We adopt as the typical example of the commonly used symmetrical 
airfoils the NACA symmetrical airfoil (ref. l4) 


T = e |i.4845^ - O. 630OX - 1.7580x^ + 1.4215x3 - 0.5075x^| 

The maximum thickness is located at x = O.3, the leading- edge radixis is 
l.le^, and the trailing- edge slope -(dT/dx) is 1.17e. The pressure 

distribution for the case e = 0.1 is shown in figxire 1. The minimum 
pressure is located at x = 0.1, and the laminar separation point, 
determined by the approximate method due to the author (refs. 15 and 16) , 
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at X = 0.61. If the transition point of the boundary layer would not 
move upstream beyond the laminar separation point, we mi^t expect to 
maintain a laminar boundary layer for more than half the surface of the 
airfoil. The flight experiments hitherto made, however, appear to give 
negative evidence for such a conjecture. 

5. Now, in order to shift the minimum pressure backward, it is 
required. to shift the position of maximum thickness (x = m) backward. 
For designing such airfoils, we represent the shapes of parts before 
and after the maximum thickness by two algebraic expressions,' For the 

forward half (O ^ x = m) 


T = e 


\j2bx + hj^x + hgX^ 


while for the rear half (m ^ x ^ l) 

T = e -^.01 + di(l - x) + d2(l - x)^ + d^Cl - x)^ 


where 


hi = 


dg = 


2 - 3^hm 
2m 

1.47 - 2di(l - m) 

(1 - m)^ 


h2 = 


d3 = 


Viim - 1 

2m^ 

dj^(l - m) - 0.98 

(1 - 


and we assign arbitrary values for three parameters, m, h (= leading- 

edge radius •*-. e^), and di (= trailing-edge slope f e) . Although the 

method has the drawback that the two expressions give different values 

of d%/dx^ at X = m, where dl/dx becomes zero, we nevertheless 
adopt it because we are in a position to vary the forward rearward 
parts most sinqily and independently. 

6. First, we fix the forward half with m = 0.5 and h = O.5, and 
vary the rear half by giving d^ the values 1.7, 2.0, 2.5, and 3.0, 

respectively. The shape of the section and the pressure distribution 
for e = 0.1 are shown in figure 2. We find from this result that, as 
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d]_ increases^ the minimum preseure point moves backward and the gradient 

of pressure rise following the minimvim pressure steepens. We also find 
that the pressure distribution in the neighborhood of the minimum pres- 
sure exhibits a wavy indentation when the value of d^ is too small or 

too large, and that there exists a certain value of d^ for which the 
pressure distribution is flat and smooth. Such a value of d^ is about 
2.5 in this case. We therefore fix the rear half with di = 2.5, and 

vary the forward half by giving h the values O.35, 0.50, O.7O, and 
1.05, respectively. The shape of the section and the pressure distri- 
bution for e = 0.1 are shown in figure 3. From this comparison, we 
find that the negative pressure bunqp immediately behind the leading 
edge decreases as h decreases, and that the maximum permissible value 
of h is about 0.7* 

The effect of thickness is shown in figure 4, in which curves of 
pressure distribution are given for different values of e, O.O6, 0.10, 
and 0.l4, but for a fixed set of parameters, m = 0.5# h = 0,5, 
di = 2.5. It is seen that the characteristics of the pressure distri- 
bution do not materially change with thickness. There is, however, a 
slight change in the pressure distribution, the maximum permissible 
value for h slightly increasing as the thickness increases. 

To see the effect of the position of maximum thickness,' we give m 
values ranging from O.35 to 0.60, varying at the same time values of d^ 
and h so that the pressure distribution becomes flat and smooth. The 
result of calculation is given in figure 5, which shows a considerable 
change in the position of minimum pressure. The change is not purely 
due to the effect of m, but it is at any rate to be noticed that the 
value of m less than 0.4 is not sufficient for shifting backward the 
minimum preseure, while increasing the value of m beyond 0. 5 is of no 
advantage, since the backward shift is then almost saturated, only the 
adverse pressure gradient being increased. 

7. From the results of calculation, we thus arrive at the conclusion 
that m must be between 0.4 and O.5 and h must be lees than O.7 in order 
that the minimum pressure occurs well downstream. Smaller values of h 
are desirable, but, on the other hand, we should like to make h as 
large as possible, because a large value of h will be advantageous in ' 
increasing the maximum lift coefficient and in preventing the inception 
of adverse pressiire gradient when the angle of attack is slightly changed. 
Even if we give h the maximum permissible value O.7, the leading-edge 
radius amounts to only 60 percent of that for the conventional NACA sym- 
metrical airfoil of the same thickness. In order to increase the leading- 
edge radius, it is required to increase the thickness, which in turn is 
accompanied by an increase in adverse pressure gradient following the 
minimum press\ire. The adverse pressure gradient should be kept within a 
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certain limit, so that it becomes necessary to make a compromise between 
conflicting requirements. Thus, we are no longer in a position to 
require the farthest possible rearward location of the minimvun pres- 
sure. We should also use a value of d^ which is somewhat smaller 

than that mentioned previously. 

Taking these requirements into acco\mt, we finally arrive at the 
design of a series of symmetrical airfoil sections, the parameters of 
which are given in the following table: 


Section 

m 

h 


Position of 
minimum pressxore 

I 

0.500 

0.35 

2.384 

0.63 

J 

.500 

.54 

1.800 

.55 

. K 

.475 

.56 

1.575 

.51 

L 

.450 

.58 

1.400 

.47 

M 

.400 

.62 

1.150 

.37 

N 

.350 

.66 

1.000 

.24 


Although section I is the most ideal for delaying the transition, in 
practice, its extraordinarily sharp nose and bl\mt tail are drawbacks. 

On the other hand, section N is too much compromised. Sections K or L 
seem to be suitable as laminar-flow airfoils for practical use. The 
ordinates of these six sections are given in table 1, while the auxiliary 
functions Bg and associated with the presavtre distribution (see 

paragraph 3) are given in tables 2 and 3> respectively. The shapes of 
the airfoil sections and the pressure distribution for e = 0.1 are 
shown in figure 6. 


DESIGN OF MEAN CAMBER LINE SUITABLE FOR LAMINAR-FLOW AIRFOIIS 


8. A symmetrical airfoil set at zero angle of attack has no lift. 

In order to obtain lift, the center line of the symmetrical airfoil must 
be curved with a suitable camber. Since the effects of thickness and 
. camber are nearly additive with regard to the pressvire distribution, the 
mean camber line which maintains the nature of the pressure distribution 
of the symmetrical airfoil will be such that it shall give a \miform dis- 
tribution of pressure difference when the thickness is removed. Evidently, 
the center of pressiire is then located at x = 0.5# so that such a camber 
line has the drawback that the travel of center of pressure is consider- 
able. To reduce the travel of center of pressure, the uniformity of pres- 
siu:e difference should be satisfied only in the foirward part of the chord. 
From the standpoint of designing the leiminar-flow airfoil, however, it is 
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required that the distribution of pressure difference is unifom 
from the leading edge to that point corresponding to the minimum pres- 
sure of the symmetrical airfoil. 

9* When the angle of attack a is small^ the expression for pres- 
sure distribution given in paragraph 3 niay be put into the form 


P 

q 


1 - 






1 + (fAg i 


Since the effect of the term (fAg + eB^) is very small, the quantity 


fG = fAc + 



a 


is required to be constant in order that the camber line shall not change 
the nature of the pressure distribution of the symmetrical section. The 
range of constancy is at least up to the position of minimum pressure of 
the symmetrical section. Putting cos I = u 2x - 1, and considering 
for simplicity the case when the minimum pressure is located at u = 0 
we prescribe that ^ 


G = constant = Gq 

for 

Vll 

0 

X ^ 0.5, -1 ^ u $ 1 

G = Gq( 1 - u^)™, m ^ 0 

for 

0.5 

^x^l, O^uSi 

figlire 7* Moreover, since 

^ = -2 nan 

1 - cos n^ 

1 

+ — 

- COS 1 

a 

1 

sin i 


sin 1 


a cannot be arbitrary, but must be so chosen that the right side of the 
equation does not become infinite at the leading edge, | = jr. It is 
given by 
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where ( *) denotes that only odd integers shoiild he 'taken for n. This 
is the so-called ideal angle of attack due to Theodorsen (ref. I?). 
Using the assun^jtion of the thin wing theory, we neglect the terms eBg 

and (fAg + We then have 


Cl = 2(1 + o^)fGo 


m 



sin2®+l| dl 


*^m0 = ■ 


1 + (1 + m)Ojn 

— Cl 

4(1 + m)(l + Ojn) 




where the moment coefficient ahout the leading edge (positive 

when nose up) at Cl = 0. Although m = 0 corresponds to the case of 

making G uniform up to the trailing edge, it seems to he in^jossihle 
to realize a finite pressure difference at the trailing edge. Moreover, 
the quantity -C^o/Cj^ (which represents the degree of center of pres- 
sure travel) is as large as 0.25 in this case. If m > 0, G vanishes 
at the trailing edge, and decreases as m increases, tending 

to 0 as m approaches 00, m = oo corresponds to the case when G = 0 
in the rear half of the chord. Increasing the value of m, however, 
steepens the press\are gradient, so the value of m from 3 to 5 seems "to 
he adeqviate. 
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Now, since 


fG sin 1 = 2 ^ nSj^ cos nl - (2 Y~ nsQ + ~ ] cos I + — 
1 \ 1 2 -kJ 2 n 


the slope of the camber line having the prescribed distribution of G 
is given by 


dM : r — 

— = 2 2_ na^ 
dx 1 sin I 


sin n| 1 


1 “ 

— 2 2 najj cos n|' 

It Jo 1 


dl’ 


COB i - cos i * 


2it 


r 


1 + 


1 + a. 


m 


log * (1 . uS)" log IjlJ!; * 

1 + U 


u 


pi (1 - u2)“ - (1 - v2)“ 

/ dv 

Jo 


U - V 




na^ 


The ordinate of the camber line is obtained hy the integration 


M= 

Jo ^ 


dx 


^ tiajj may be determined by the condition that M = 0 at x = 1 . We 
1 

call Dm the camber line thus determined. The equations for camber 
lines for m = 0 , 1, 3, 5^ and <», namely Dq, Di, D3, D5, and D„, 
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are given below, their in^ortant characteristics being summarized in the 
table. It is to be noted that f is the maximvun value of M, e is 
the absolute value of the zero-lift angle, and e and a are measured 
in radians . 


m 

o./Cl 

€/Cl 

f/CL 

-CffiO/CL 

0 

0 

0.1592 

0.0552 

0.2500 

1 

.0380 

.1211 

.0711 

.1750 

3 

.0609 

.0983 

.0790 

.1213 

5 

.0703 

.0888 

.0816 

.0979 

CO 

.1103 

.0488 

.0874 

0 



1 

2 log 2 


1^(1 - u) log (1 - u) '+ (1 + u) 


log (1 + u) 


} 


i (5 + u) log 2 - — (1 + u) log (1 + u) + ^ u3 log |u I - 
Cl 5 5 5 ' ' 

i (1 - u)^(2 + u) log (1 - u) + i (1 - u^) 

5 ' , .5 

■ ^ ^3(35 T 21u^ +: 5u^y log . |u ^ 

~ (1 r. u)^(i6 + 29u + 20u2 + 5u3) log (l - u) + 

51 , • 

: i:(l - I^)(i 76 - 8lu - .l-^u^ + 30 u3 + 60 u*^) 
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® 5 * ^ (9*^9 + hsin) log 2 - ^ (1 + u) log (1 + u) + 

Cl 9^9 949 

^ u3(ii55 - 1386u2 + 990u'^ - 385u^ + 63 U®) log |u| - 

^ (1 - u)6(256 + 843u + 12i8u 2 + 938u3 + 3 y8u^ + 

63 u^) log (1 - u) 4 

31680u3 4 68792u^ - 
D„: I = |(1 4 u) log 2 - 

^ log 3 I 


— (1 - u2) (35072 - 28535U - 66088 u2 + 
113880 ' 

17430u5 - 36120u^ 4 3780u'?' 4 7560u3) 


(1 4 u) log (1 4 u) 4 U log ju 


Shapes of these camher lines are shown in figure 7, the ordinates 
of them are given in table 4, and the auxiliary functions and Ag 

(see paragraph 3) and the pressure difference distribution G are given 
in tables 6, and 7> respectively, 

10. The calculation made previo\isly is only approximate, neglecting 

the thickness. It is therefore desirable to check the result by actually 
calculating the pressure distribution for the specified angle of attack 
taking both camber and thickness into account. As an example, we con- 
struct an airfoil by applying the thickness form K with e = 0.15 nor- 
mal to the chord around the camber line with f = 0.02 (the 

resulting airfoil is designated as - 2OI5). We calculate the pres- 

sure distribution by the formula of paragraph 3 for the optimimi design 
condition a .= 0.99° and Cl = 0.245. The result is shown in figure 8. 

The nature of the pressure distribution remains similar to that of the 
symmetrical airfoil, so we may consider that the approximate determina- 
tion neglecting thickness gives results sufficiently accurate for practi- 
cal purposes. 

11. In designing the camber line Dj^, we have assumed for simplicity 

that the pressure difference G is constant for u ^ 0. This corresponds 
to the case when the symmetrical airfoil has its minimum pressure in the 
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neighborhood of u = 0. Therefore, the camber line Djjj is adequate to 
be combined with the s^nmetrical section J or K. If, however, the 
symmetrical section is adopted in which the minimum pressure is located 
finrther upstream, it is not only not necessary to maintain G constant 
up to u = 0, but also of disadvantage because it makes it difficult to 
reduce the value of -Cmo/CL. 


To reduce the range over which G should be maintained constant, 
we may proceed in the following way. Assuming for instance that G 

should be constant from u = -1 to u = - and using a new variable 
uq = ^ (1 + 3u) , we prescribe that 


G = Gq (1 - ui2)3 for - i S u < 1, 0 < ui < 1 

The calculation may be performed similarly to the case of Djj. The 
resulting camber line is designated as F^. The camber line lying in 
the middle between and is also designed, and designated as Ejj^, 

Their important characteristics are given in the following table together 
with those of D^. Other numerical data for these camber lines are given 
in tables 4 to 7* 


Camber line 

a/CL 

e/CL 

f/CL 


% 

^3 

0.0703 

.0752 

.0764 

0.0888 

.0840 

.0827 

O.O816 

.0813 

.0795 



Since the camber lines Ei|. and F^ enable us to maintain G con- 
stant up to the point x = 0.42 and x = O.33, respectively, they are 
adequate to be combined with the symmetrical sections L and M, respec- 
tively. The pressure distribution is shown in figure 10 for the airfoil 
obtained by applying the thickness form ,M with e = O.15 around the 
camber line F^ with f = 0.02. The optimum design condition corresponds 

to a = 1.10° and Cj, = 0.252. 
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EXPERDffiNTS ON LAMINAR -FLOW AIRFOIIS 


12. In order to ascertain whether it is possible to prevent the 
forward movement of the' boundary layer transition by shifting the mini- 
mum pressiire on the airfoil surface, we have to perform experiments in 
a low turbulence wind tunnel or on the actual airplane in flight. When 
the Reynolds number is not too large, however, we can still use a con- 
ventional wind tunnel in which the stream turbulence is relatively small. 
So we made at first conqjarative measurements on two symmetrical airfoils, 
NACA 0010 and L.B, 2 h in the I.5 m wind tunnel of the Aeronautical 
Research Institute. L.B. 2 k is a laminar-flow airfoil of 10 percent 
thickness, already shown in figure 3. The theoretical pressure distri- 
bution is also given in figure 11. The minimum pressure is located at 
X = 0.64, and the laminar separation at x = O.77. The wind tunnel was 
of the lowest turbulence level available for the author, the critical 

Reynolds number of a sphere being 3.66 x 10^ and the transition Reynolds 

number of a flat plate I.05 X lo3 (see paragraph 28). In order to raise 
the Reynolds number as high as possible, unusually large models were 
used. They were made of laminated mahogany, of highly polished surface, 
of 0.8 m span, of 1.2 m chord, and fitted with end plates I.3 m X 0.6 m. 
Since the model was large con5)ared with the size of the tunnel and the 
end plates were not sufficiently large, the results for a given airfoil 
may not correspond even approximately with those for the same airfoil in 
an undisturbed two-dimensional flow. Our object, however, was merely 
to ascertain the relation between pressiire distribution and transition, 
and it seemed reasonable to expect that the relation will not be seriously 
affected by limitations in the conditions of the experiments. As a matter 
of fact, marked difference was found in the calculated and measured dis- 
tributions of pressure, the latter of which was measured along the median 
section of the model with a static tube of 1 mm diameter (fig. 12). 2 
This discrepancy, however, is immaterial, since our object was merely to 
con^iare the two airfoils, both of which are affected quite . similarly by 
experimental limitations. 

13. The angle of attack of the model was zero, and the wind speed 
was varied from 6 to 40 m/s. The local drag of the median section was 
determined from wake measurements, that were made in the section 11 cm 
behind the trailing . edge. Measurements of static and total pressures in 
the wake were made, respectively, with a static tube of 2.5 mm external 
diameter and a pitot tube with a flattened mouth of O.65 mm external 
depth and 2.6 mm width. The profile drag coefficient CpQ was obtained 


2 

The measured values are those for a Reynolds number of about 2 X 10^. 
The distribution of pressure changes but little with the Reynolds number. 
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from the measured pressures hy Jones* formula (ref. I 8 ) . Figure I 3 ^ 
presents CpQ plotted against Reynolds nvuober R referred to chord 

length. For a lower range of R, the drag of L.B. 24 is higher than 
that of NACA 0010, the reason prohahly being that a tiorhulent boundary 
layer associated with the laminar separation is established at a higher 
Reynolds number for the former airfoil than for the latter. For a higher 
range of R, however, the condition is reversed, L.B. 24 giving a drag 

less than half that of NACA 0010 for R higher than 2 X 10^. This is 
probably due to the fact that the transition may occvir much later for 
L.B. 24 than it does for NACA 0010, as also obseirved from the con^iari- 

spn of wake conditions for the two airfoils (fig. l4).^ 

14. In order to verify the aforementioned supposition, a pitot tube 
with a flattened mouth of external depth O .9 mm and width 2.7 mm was 
placed in contact with the airfoil surface, and the wind speed, and con- 
sequently the Reynolds number R, were determined at which the indicated 
total pressure G* divided by the dynamical pressure q of the undis- 
turbed stream begins to rise suddenly. The results are shown in fig- 
ure 15 . From this figure, the dependence of the transition point on 
Reynolds nvunber as shown in figure I 6 is obtained. At the same Reynolds 
number, the transition occurs much farther from the leading edge for 
L.B. 24 than for NACA 0010. Even at the highest Reynolds number reached, 
L.B. 24 has a transition as far back as x = O. 8 O. This is somewhat 
beyond the laminar separation point, x = 0.77# which is calculated from 
the theoretical pressure distribution. However, this is not contradictory, 
because the actual pressure distribution differs from the theoretical one 
in a manner to delay the transition (fig. 12). 

15 . With further increase in Reynolds number, the transition may 
move toward the leading edge, but it seems improbable that the transition 
moves forward beyond the minimum pressure. It is highly desirable to 
check this point also by wind tunnel experiments, but all the wind tunnels 
now available to the author are of no use for' making measurements at 
sufficiently high Reynolds numbers, because the transition is prematurely 


^In this figure, the- curves L and T represent the drag of a 
flat plate when the boimdary layer is entirely laminar and entirely 
turbulent, respectively. The curves NV and NF represent the drag 
of airfoil NACA OOO 9 measured in the NACA Variable-Density Wind Tunnel 
and NACA Full-Scale Wind Tunnel, respectively. 

^In this figure G and p are the total and static pressures in 
the wake, respectively, and Gq is the total pressure outside the wake 
all being measiured from the static pressure of the undisturbed stream, 
y is the distance across the wake, and t is the chord of the model. 
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induced by the turbulence of the stream (see paragraphs 27 to 29 ). It 
seems urgent to build a special wind tunnel of low turbulence level 

simplest to rely upon experiments i^ 
actual Such a hope of the author was fortunately realized by 

the specially planned flight ejq)eriment, which was performed at the 
Navy Aeronautical Technical Arsenal (ref. I 9 ) . 

^■p +>, 0^1 airplane used for the experiment was a biplane; two portions 

+? covered with the airfoil to 

be tested. The test portions were of chord 2.4 m, made of Japanese Hinoki 

fitted with a partition fence of small height at " 
bn+?n« ^Z° portions were placed symmetrically, pressure distri- 

wh^fthf a perfomed on the starboard portion, 

while the boundary layer was observed on the port portion. The airfoil 

appropriate design now considered, because 
t was required to put it on the original section of the airplane, and 
TOreover, to determine the section before completion of the final desim 
calculation. It has the following characteristics: 

Mean camber line: M = 0.0667x(l - x)(l - x + x2), f = 0.0125 


Thickness distribution: 


e = 0.12, m = 0 . 45 , h = O. 56 , d.x = I .60 


The camber line is similar to Dq of paragraph 9, but there exists a 

unifomity of G in the neighborhood of leading and 
trailing edps. The thickness distribution is similar to L of para- 
graph 7 , but the trailing edge slope is somewhat larger than L. 

17« Results of flight experiments are summarized in figure I 7 
and 18. In figure I 7 , the section lift coefficient Cl, obtained by 
inte^ating the pressure distribution curve, is shown by a broken line 
plotted against the Reynolds number R referred to the flight speed 
and chord length, and Cl is again shown by a solid line plotted against 

the profile drag coefficient C^q determined from the wake measurements. 

the, measured pressure distribution is shown in comparison 
vitn the theoretical one (two-dimensional potential flow) having the 
same value of Cl. The transition points estimated from the change in 
boundary layer velocity profiles are also marked. Generally speaking 
the measured pressure distribution agrees fairly well with the theo- 
retical one, although a slight difference appears when Cl becomes 
large. An adverse pressure gradient is found on the lower surface when 
L 8 small, thus resulting in the transition point being observed 
unexpectedly far foiVard. Such a discrepancy in pressure distribution 
as conqpared with the theoretical one seems to be probably due to the 
fact t^t the span of the test portion was not sufficiently large. As 
a res , e profile drag coefficient Cj)^ has the minimum value 0.0042 
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at about Cl = 0.26, which is larger than the value Cl = O.I8 theo- 
retically estimated on the assumption that the transition occurs far 
back on both upper and lower surfaces. Therefore, the observed value 
of Clq, although much smaller than that of the conventional airfoils, 

seems to be still somewhat large when compared with the optimum case. 

At any rate, however, no transition was found to occur upstream of the 
minimum pressure. It is important to note that such experimental evi- 
dence was obtained on an airfoil section in which the minimimi pressure 
is located further downstream than on the conventional one . This finding 
will give valuable data to establish a basis for design of the laminar- 
flow airfoils. 


ESTIMATION OF THE DRAG OF LAMINAR-FLOW AIRFOIIS 


18. As mentioned previo\isly, the results of flight experiments 
seem to support the basis for the design of laminar-flow airfoils, 
namely, the possibility of maintaining the boundary layer laminar at 
least up to the minimum pressure point. It is interesting, therefore, 
to estimate the drag of laminar-flow airfoils by assuming a laminar, 
boundary layer from the leading edge to the minimum pressure point and 
a turbulent' boundary layer downstream to the trailing edge. 

For the laminar boimdary layer, the momentum thickness is given by 


02 


O.kkv 





ds 


with a sufficient approximation (ref. 16) , where u^ is the velocity 

outside the boundary layer and s is the distance measured along the 
airfoil surface from the forward stagnation point. Writing t for the 
chord length and V for the velocity of the vindisturbed flow (velocity 
of flight) , and putting 


u-j^ = UV 


s = at 


R = — 


we have the nondimensional expression 
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£sl = 2iii r= 

^2 RUa Jo W 


where the subscript a refers to the point of minimum pressure. 
Applying then the solution due to Buri (ref. 20) for a turbulent 
boundary layer assimed to extend from the mlniTnnm pressvire point to 
the trailing edge, we have the resvilt 


3 Q 


where the subscript b refers to the trailing edge. The numerical 
values originally given by Buri are slightly modified so as to agree 
with measurements when applied to the flat plate. 

According to Squire and Young (ref. 21), the profile drag coeffi- 
cient is given by 


% = I («U - 

Where the subscripts u and Z refer to the upper and lower surfaces, 
respectively. The exponent 3«2 of Ujj has been obtained by assuming 

the ratio of displacement and momentum thicknesses equal to 1.4, But the 
ratio seems to exceed 1.4 near the trailing edge, so we replace 3.2 by 
3.4 with a view to ii^roving the accuracy and at the same time to 
sia^Jlifying the algebra. Since 
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we have 



O.OTl^R 


- 1/5 




(ljR-3/8 + t^^V5 


where 





If the velocity distrihution = UV is calculated hy assuming the 

potential flow of an ideal fluid, it is desirable to modify the distri- 
bution to take account of the effect of separation near the trailing 
edge. We tentatively modified the distribution of U such that the 

Burl parameter T = (e/uj) (dui/ds) at the trailing edge for 

the case when the boundary layer is assumed turbulent from the leading 
edge, namely 


U^6 \da Jo 

shall not become smaller than -O.06. In almost all the cases, values 
of Uij thus modified are found in the range between O.95 and 1 . 00 . 

19. Applying this method of calculation, the profile drag coeffi- 
cient is estimated first for a series of symmetrical airfoils set 

at zero angle of attack. The, series consists of the six symmetrical 
airfoils, I, J, K, L, M, N, as given in paragraph 7 and the NACA con- 
ventional airfoil. Values of Cd^ at R = 2 X 10 for three different 

thicknesses (maximum thickness in terms of chord e = 0 . 10 , O.15, 0 . 20 ) 
are shown in figure I9 plotted against the position of minimum pressure, 
seems to decrease almost linearly as the minimum pressure is shifted 

backward, the most ideal airfoil I giving a value about half of that of 
the NACA conventional airfoil. If it is desired to realize a profile 
drag of two-thirds of the conventional airfoil, it will be required to 
use the symmetrical airfoil L with the maximiam thickness located at 
45 percent chord from the leading edge. 
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Effect of camber is relatively small. If, for instance, the 
center line of the symmetrical airfoil K vith e = O.I5 is curved 
into the camber line D5 with f = 0.02 (see paragraph 9), the esti- 
mated increase in profile drag at the optimiim angle of attack is only 
0.0001. For f = O.OU, it is O.OOO3. 

Finally, we compare the laminar-flow airfoil with the most exten- 
sively used airfoil, NACA 23012, for which the leading-edge radius is 
0.0158, and the optimum lift coefficient corresponding to the minimum 
profile drag coefficient is about O.I5. If we consider the symmetrical 
airfoil section K (h = O.56) combined with the camber line D^, it is 

necessary to use the thickness e = O.I5 in order to obtain the same 
magnitude of leading-edge radiiis, and the camber f = 0.012 in order to 
realize the optimum lift coefficient 0.15.^ Therefore we construct an 
airfoil by applying the symmetrical form K with e = 0.I5 normal to 
the chord around the camber line D5 with f = 0.012. We call it 

- 1215. The angle of attack corresponding to Cl = O.15 is 1.56° 
for NACA 23012 and 0.60° for - 1215. The pressure distribution for 

that condition is shown in figure 20. 

We then estimate Cj)q for the two airfoils by the method explained 

previoiasly. The results are shown by broken lines in figure 21. In 
order to check the results, measured values taken from various sources 
for the two airfoils and similar airfoils are also plotted in the same 
figure by different marks. The mark o refers to the value obtained 
by flight experiments on a smooth surface, and • refers to that 
obtained by wind tunnel experiments where the stream turbulence has no 
effect on transition. The mark + refers to the flight experiment on 
a rough surface, while x refers to the wind tunnel experiment where 
the stream turbulence causes the transition to occur prematurely. There- 
fore, only o and • are adequate for our present pxirpose. Drawing 
curves^ through these points and extrapolating to higher Reynolds numbers, 
we find that the result agrees fairly well with the estimated values. 
Therefore, we may consider that the method of estimating Cp^ ie suffi- 
ciently accurate at the Reynolds numbers corresponding to actual flight 


5lhe calculation developed in paragraphs 8 to 11 refers to the poten- 
tial flow of an ideal fluid, so that it gives the slope of lift curve 
dCL 

— — = 2if. In real fluids, however, the slope of lift ciorve amounts to 
da 

only 80 to 90 percent of the theoretical value. If we take this effect 
into account, we have to increase the necessary amount of f by 10 to 
20 percent in order to realize the given lift coefficient. However, 
such a slight change in the value of f will scarcely affect the esti- 
mation of Cpq . 
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conditions. Con5>ari8on of two airfoils, laminar -flow and conventional, 
also suggests the possibility of ^iO percent reduction in profile drag 
by using a fairly practical laminar-flow airfoil. 


CONSIDERATION OF THE AIRFOIL WITH UNIFORM DISTRIBUTION OF PRESSURE 


20. The fact that, so far as flight experiments with smooth wings 
are concerned, the boundary layer transition occurs only in the region 
of rising pressure, not only warrants the principle- of designing the 
laminar-flow airfoil by shifting the minimum pressvire backward, but 
also suggests the possibility of delaying the transition by using an 
airfoil with tmiform distribution of pressure. Therefore, in para- 
graphs 21 to 23 , the shape- of such a symmetrical airfoil is determined 
by a method similar to that used for designing the camber line of inmiTinT 
flow airfoils, and the airfoil was examined by wind tunnel experiments. 

In paragraphs 24 to 26, a calculation is made to inquire about the method 
of sucking away the boundary layer over the region of rising pressure in 
such a way that the boundary layer velocity profile shall remain the same 
as that for the point of minimum pressure. 

21. Consider the syimnetrical airfoil set at zero angle of attack. 
According to the formula of paragraph 3> the pressure distribution is 
given by 



where p is the pressure acting on the airfoil evurface^ measured from 
the static pressure of the undisturbed stream, q is the dynamic pres- 
sure of the undisturbed stream, x = i (l + cos I) is the coordinate 

along the chord, and the half -thickness of the airfoil is ejqjressed in 
the form 


00 


1 


bjj sin n| 
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If the thickness is sufficiently small, the square of 


dx 


00 

= “2 H 


1 


cos n| 
sin 4 


may he neglected, so that the condition of uniform distribution of pres- 
sure is satisfied by putting all the coefficients b^, other than b]_, 

equal to zero, namely, by an elliptic section. In order to take the 
thickness into account approximately, we substitute the value of dT/dx 
for the elliptic section into the denominator of the expression for p/q. 
Then, writing e for the maximum thickness in terms of the chord we 
get " 


sin I + 2 ^ nb„ sin n| 



sin^l + e^cos^l 


hence 


2 nbj^ sin n 4 = - sin I + B\|sin^l + e^cos^l 
1 

1/2 

where B is the constant value of (l - p/q) . We have therefore 
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sin i 


dx 



sin nl' 

1 


sin I' dt' 
cos I - cos 


where 




f 


_ 1 

r ■ 

sin - B^sin^^ ' 

+ e^cos^^' * 

It 

Jo 






- -- 

= u 

B 

— — ^ 

2ku cos"^e + 



Tt 




Bin I’ dl* 
cos I - cos I' 


— ^ ( 1 + u) (1 - k^u + e^ll - k^u^) 

1 - k^u^ log ^ 

(1 - u)(l + k^u + eil - k^u^) 

' -- 


u = cos I = 2x - 1 


k = 



e 


2 


Upon integrating we get 



The integral is evaluated hy a numerical method, and the value of the 
constant B determined from the condition that ^ ^ when u = 0. 

The numerical results for three values of e are given in table 8, T 
and dT/dx being expressed in terms of those for the elliptic section. 
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The shapes of the airfoils are shown in figure 22. The shape resembles 
an elliptic section, although it is somewhat fuller at the ends. It will 
also be seen that the constant B, as shown in the following table, is 

slightly smaller than 1 + e, the maximum value of (l - p/q)^^^ for the 

elliptic section. Since the values of (dT/dx) (-eu/lBT - u^) are not far 
different from 1, it seems to be sufficiently accurate to .substitute 
the value of dT/dx for the elliptic section into the denominator 
of 1 - p/q. 


e 

B 

0.1 

.2 

.3 

1.097 

1.188 

1.273 


22. The uniform distribution ^pf pressure requires, however, an 
uifinxte pressure gradient at both leading and trailing edges. In order 
to see to what degree such a sharp pressure gradient may be realized in 
actual fluids, measurements were made bn a model of the airfoil section 
with uniform distribution of pressure with e = 0.1 (we call it U.P. 0010) 
in the I.5 m wind tunnel of the Aeronautical Research Institute. The 
model was made of laminated mahogany, of 0.8 m span, of 0.8 m chord, and 
fitted with end plates I.3 m X 0.6 m. Measurements of pressure distribu- 
tion, wake traverse and boundary layer transition were similar to those 
already mentioned in paragraphs 12 to l4. 

The pressure distribution along the chord is shown in figure 23 
for three values of R, the Reynolds number referred to chord length. 

The observed value is somewhat high conqpared to the theoretical 

value £ = -0.203, the discrepancy probably being due to the excessive 

size of the model in proportion to that of the wind tunnel. At any rate, 
however, the pressure distribution is nearly uniform. The ’lack of uni- 
formity exists at both edges due to the impossibility of realizing the 
infinite pressure gradient. The boundary layer separatee near the 
trailing edge, but the effect of separation becomes small as the Reynolds 
number increases . This scale effect seems to be of the same natiu*e as 
that responsible for the sudden drop in sphere drag; the boundary layer 
separates in a laminar state when the Reynolds number is low, while it 
becomes turbulent before separation when the Reynolds number is high, 
thus being able to proceed against a larger pressure gradient. This is 
also seen from the measurements in the wake, where the indentation of the 
curve of total pressure distribution is shallow and wide for low Reynolds 
numbers, while it becomes deep and narrow as the Reynolds number increases. 
As a result, the profile drag coefficient decreases considerably 
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as the Reypolds niunber increases, as shown in figure 2k, Transition to 
turhnlence was found very near to the trailing edge, occurring downstream 
of X = 0.9 in the range of measurements. Measured values of the drag 
of the model when a piano wire of 0.5 mm diameter was placed at x = 0.8 
and X = 0.9, respectively, are also plotted in the same figm:e. The 
drop in drag occurs at a lower value of the Reynolds number when the 
svirface is rou^ened by the wire. 

23. The profile drag coefficient of the airfoil U.P. 0010 is shown 
below in comparison with other symmetrical airfoils at R = 2.2 X 10^; 


NACA 0010 
L.B. 24 
U.P. 0010 

U.P. 0010, a wire at 


X = 0.9 


= 0.0064 (fig. 13) 
0.0032 (fig. 13) 
0.0059 

0.0044 


All the airfoils are of 10 percent thickness. NACA 0010 is a conventional 
airfoil, and L.B. 24 is a laminar-flow airfoil with far back minimum 
press\ire. The drag of U.P. 0010 is between that of these two airfoils, 
the drag when a wire is placed being nearly the mean of the two. This 
resxilt seems to be interesting in that the drag of an airfoil yith a 
blunt tail is smaller than commonly considered. 

The airfoil with uniform distribution of pressure will also probably 
be favorable when used at high subsonic speeds. Even if the shock wave 
occurs at high subsonic speeds, the increase in drag will remain small 
when the boundary layer does not separate. This expectation was really 
verified by the experiment due to Kawada and Kawamura (ref. 22), the 
drag of the airfoil U.P. 0010 being smaller at high Mach numbers as 
con5)ared with other airfoils. 

24. From the fact that the boundary layer transition occurs only in 
the region of rising pressure, we may also expect suction of the boundary 
layer to delay transition. For example, if the boundary layer is sucked 
into a slot, there is a well-known sink effect (ref. 23) which relieves 
the adverse pressure gradient somewhat upstream of the slot. We may con- 
sider an alternative possibility. That is, we assume that the boundary 
layer is laminar in the region of falling pressure, and that it remains 
laminar also in the region of rising pressure provided that the boundary 
layer profile is the same as that at the minimum pressure point. We 
then ask what suction arrangement must be applied in order to realize 
such a condition. 


25. We denote by s the coordinate meas\ired along the surface, 
y perpendicular to the surface, & the boundary layer thickness, u 
the velocity in the boundary layer, u^ the velocity outside the boundary 
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dp 

ds ^ the pressure gradient, and tq = -h(Su/ay) ^ 

the skin friction* at the surface. Assuming the surface (y = O) is made 
porous, through which the fluid is sucked with the velocity c, we have 
the equation of continuity 


w - c = 


'T 

ds Jq 


(a) 


and the equation of momentum 


puj^w - p 


d r 2 j c dp 

— / u^ dy = tq + 6 — 

ds Jo ’ 


ds 


Cb) 


where w is the velocity of fluid entering the boundary layer through 
»y = 6." The equation of motion reduces to 


■ W/,.o “ 


(0) 


for y = 0. 

Now, the velocity profile in the boundary layer may be approximated 
by the Pohlhausen polynomial 





(d) 


when neither pressure gradient nor suction exists. For this profile we 
have 



(e) 


kkl 


(f) 
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To simplify the calculation, we assume that the expressions (d), (e), 
and (f) still hold when hoth pressure gradient and suction exist. Then 
we have from (h) • 


^ ^ -L 

ds 37 367 ds 


Integrating we have 


m 6^ 



+ 


3^60 P® 
37 JsQ 


ds 


(g) 


(h) 


V, 208 ^ . . 

where m - ^0 initial position of suction, which is 

the minimiom pressure point in the present case. We have also from (a) 

l84l ^^1 ^ , 

c 6 (i 

3670 ds 


If we substitute (d) and ( i) into both sides of the equation (c), however, 
the left and right sides become 1.003ui^dui/ds) and ui^dui/dsj, respec- 
tively. This contradiction is evidently due to the criode assunq)tion of 
using (f) in spite of the presence of pressure gradient and suction, but 
we may overlook the error becaiise it is small. 


26. We apply the calculation to the symmetrical laminar-flow air- 
foil of 10 percent maximimi thickness, L.B. 24, set at zero angle of 
attack. The velocity distribution u^/V calculated from the potential 

flow of ideal fluids is used,^ the maximum velocity (minimum pressure) 
being located at 64 percent of the chord from the leading edge ( sq = 0.6 

Applying a distributed suction downstream of the minimum pressure point 
so as to maintain the velocity profile in the boundary layer the same as 
at this point, we have the boundary layer thickness 8 the required 

^It is assumed that the velocity distribution is not affected by 
the suction. Theoretically u^ should be 0 at the trailing edge, but 
the distribution was somewhat modified so as to give up = O.85V there. 

The effects of these assimiptions appear to be too small to affect the 
result materially. 
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suction velocity c as shown in figure 25. It is to be noted that s 
is measvired along the surface from the leading edge, and R is the 
Reynolds number based on the chord length t and the velocity of the 
imdisturbed stream V. Integrating the area under the cvirve of c, we 

have the total amount of suction 1.3 vi/R per unit span of the two s\ir- 
faces. If we assume span = 35 m, t = 5 m, V = 200 m/s, v = 0.I5 cm^/s 
the total amount of suction amounts to 5.6 m^/s, which will require an 
exit area of only 0.028 m^ vhen discharged with the velocity eqiial to V 
Integration of tq gives the drag coefficient Cdq = O.OOO5. This value 

may be compared with Cdq = O.OOO3 for the flat plate with laminar 
boundary layer, = 0.00i^4 for the flat plate with turbulent boundary 

layer, and Cdq = 0.0025 for L.B. 2 k without suction. If the thickness 

of the airfoil is doubled (20 percent chord), then the amount of suction 
w'ill be nearly doubled; the drag is however almost unchanged. 

It should be noticed again that the calculation is based on the 
assunqjtion that no transition occurs if the velocity profile in the 
boundary layer maintains the form at the minimum pressure. It is the 
purpose of the calculation to show that extraordinarily low profile drag 
may be expected with a relatively small amount of suction under such a 
condition. 


PREMATURE TRANSITION OF BOUNDARY LAYER - EFFECT OF STREAM TURBULENCE 


27. Although the transition of the boundary layer occurs only down- 
stream of the minimum pressure , so far as flight experiments on smooth 
airfoil siufaces are concerned, there are many exan^iles of wind tunnel 
experiments in which the transition moves upstream of the minimum pres- 
sure. This seems to be due to the prematxare transition caused by the 
turbulence in the wind tunnel stream. For exan5)le, the transition on 
the airfoil L.B, 24 was found only downstream of the minimum pressure 
in the range of Reynolds numbers covered by the author ’s wind tunnel 

experiments (the Reynolds number based on chord length up to 3 ^ 10^; 
see paragraphs 12 to l4) ; as a result very low values of the profile drag 
coefficient Cpq were observed. The same airfoil, however, when tested 

with a larger model of 2 m chord in the 2. 5 m wind tunnel of Kawasaki 
Aircraft Company, Glhu, gave the result as shown in figure 26, in which 

CpQ increases considerably when the Reynolds number exceeds 5 x 106 

(refs. 24 and 25). There is reason to believe that the increase in drag 
is due to the effect of stream turbulence. The boundary layer observa- 
tion at the Kawasaki wind tunnel shows that the transition is found at 

50 percent chord (x = 0.5) for the Reynolds number 6 x 106 and moves 



28 


NACA TM 1351 


further forward as the Reynolds number increases. The bo\indary layer 


velocity profile observed at transition has a form factor 


6f; 

0 


2.6 


(8* and 0 are the displacement and momentum thicknesses of the 
bovindary layer ) , which is very near to the value for the case of zero 
pressure gradient. This result seems to suggest that the transition is 
at least not correlated with the laminar separation (ref. l) . 


28. In order to verify this conjecture, it is desirable to show 
that the transition in the boxmdary layer along a flat plate occ\irs 
under the same condition, because that transition may be considered to 
be independent of the laminar separation. Unfort\jnately, however, no 
flat plate was measured in the Kawasaki wind tunnel. Therefore we pro- 
ceed in a somewhat indirect way. We assume that the degree of stream 
turbialence is represented by the conventional critical Reynolds number 
of the sphere, Rq, and the condition of transition due to turbulence 

represented by the local Reynolds number, Rq = -i— , at transition on a 

V 

flat plate, where is the velocity outside the boundary layer' and 

0 is the momentum thickness of the boundary layer. It is generally 
accepted that the turbulence in the wind tunnel stream will give a 
fluctuation of pressure gradient, as a result of which an instantaneous 
and intermittent separation will occur. Such an instantaneous and inter- 
mittent separation, however, does not necessarily lead to the transition 
into turbulence; for the transition really to occur, it seems probably 
necessary that the Reynolds number Rg which represents the ratio of 

inertia pu^^^ to viscous stress pva-^/Q exceed a certain critical value. 

It is also expected that the critical value depends on the degree of tur- 
bxilence; it must increase as Rq increases. This is really shown by 
the experimental data hitherto published, which are given in the following 
. table and also by white circles (o) in figure 27. The available data are 
scanty, especially because the experiment on a flat plate is very diffi- 
cult. It was necessary for the author to perform a new experiment (ref. 30) 
with a view to adding one point in the range of high Rq. 


— 

Rc 

R0 

Wind Timnel 

Reference 

l.Uo X io5 
2.75 X io5 
2.20 X io5 
3.66 X io5 

0.21 X io3 
.70 X io3 
A2 X io3 
I 1.05 X io3 

1 

National Bureau of Standards 
National Bureau of Standards 
N.P.L. Conqjressed Air Tunnel 
Aero. Res. Inst. I.5 m Tunnel 

26 

26 

27, 28 
29, 30 
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29. Now, we calculate the value of Rg at transition of the air- 
foil L.B, tested in the Kawasaki wind tunnel, and correlate it with 
the critical Reynolds number Rc of that tunnel. We analyze similarly 

the other available data, and summarize the result in the following 
table. The values of the form factor 6 */ 0 , not shown, were all found 
in the range from 2.1 to 2.7. We then plot the data also in figure 27 
by black circles (•) . The black circles are seen to define a single 
curve together with the white circles already mentioned. This result 
seems to support the supposition tliat the transition under consideration 
is mainly caused by the stream turbulence, but not correlated with the 
laminar separation. 



% 

Model 

■ 

Wind Tunnel 

Reference 

1.85 X io5 

2.10 X 105 

3.50 X io5 
3.65 X io5 

O.U X 103 
.56 X io3 
.95 X io3 
1.08 X io3 

Symmetrical airfoil 
Airship model 
Airfoil N-22 
Airfoil L.B. 2 k 

N.P.L. 7 ft 
M.I.T. t| ft 
NACA Full-Scale 
Kawasaki 2.5 m 

1 

31, 32 

33, 34 
12, 35 
24, 25 


In reference 12 (the third line in the preceding table), the same 
airfoil was examined both by the full-scale wind tunnel and by the flight 
test^. We calculate the form factor b */6 from these tests and plot the 
values against s/t in figure 28, where t , is the chord length and s 
is the length measured along the surface from the leading edge. The 
value of 6*/e at transition is 2.6 in the wind tunnel, while it reaches 
as high as 3.1 and drops sharply in the flight test. The minimum pres- 

s 

sure is located at ^ = 0«l8, and the laminar separation calculated from’ 

the measured distribution of pressure at = 0 . 3 ^» This example is veiry 

interesting because the cause of transition is quite different in the tvo 
cases (namely, it is due to the stream turbulence in the vind tunnel, 
while it is related to laminar separation in the flight test), although 
the positions of transition are almost the same. 


PREMATURE TRANSITION OF BOUNDARY LAYER - EFFECT OF SURFACE ROUGHNESS 


30. Up to this point, we have only considered the case when the sur- 
face of the airfoil is smooth. If the surface is rough, however, there 
is a possibility that the transition may also be caused prematurely by 
surface roughness. So, it is in5)ortant in practice to estimate the 
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approximate order of magnitude of the permissible roughness in the lami- , 
nar boundary layer. Nothing has been known concerning this problem^ 
except a mere conjecture or fragmentary data. Schiller (ref. 36) sug- 
gested that a local separation occurs and hence leads to transition 
when the Reynolds number kuj^V exceeds a certain critical value ^^crit^ 

where k is the height of projection and u^ is the velocity at the 

top of projection. The exact value of Rcrit known^ but it is 

not likely to differ much from the critical value of the Reynolds number^ 
above which vortices are shed from the obstacle of the same shape as that 
of the projection placed in a uniform stream. The experimental result 
due to Wieselsberger (ref. 37) shows that such a critical Reynolds number 
is roughly 50 for a circular cylinder. Assuming that the height of pro- 
jection k is small, and that the presence of the projection in no way 
alters the character of the flow, we have the shearing stress at the 

surface Tq = pv(u}^/k) . Using the so-called friction velocity v^ = iFo/p 


instead of u^, we have then permissible roughness 

"is therefore given by kv^/v = or, with R^rit =50, kv^/V = 7. 

On the other hand, according to Nikuradse’s experiments on roughened 

kv 

pipes (ref. 38), the critical Reynolds number is — ^ = 4, above which 


the roughness projections disturb the laminar sublayer of the turbulent 
boundary layer and hence increase the pressure drop. It appears there- 
fore that the permissible roughness, is smaller in the turbulent boundary 
layer as compared with the laminar boundary layer. This is confined by 
a British flight experiment (ref. 39) on the airfoil section of 10 feet 


chord (Reynolds number 1.8 X lOT) ^ because the effect of camouflage paint 
of 0.001 inch thickness increased the drag by about 6 percent without 
moving the transition forward. At any rate, however, such a estimate is 
nothing but mere conjecture. With a view to making the estimate more 
definite, we performed wind tunnel experiments, although of small scale, 
of quantitative character (paragraphs 31 to 32). 


31. A polished aluminum plate, 80 cm long, 60 cm wide, and 3 nmi thick, 
was held horizontally in the 1.5 m wind tunnel of the Aeronautical Research 
.Institute. So that the flow at entry would not be disturbed, the leading 
edge of the plate was rounded, and the plate slightly tilted so that the 
forward stagnation point was on the same surface as that where the observa- 
tion was made. The tilting, however, was so slight that the static pres- 
sure was observed to be practically uniform along the plate. The plate 
was roughened by a wire, which was stretched across the flow, in contact 
with the plate. The diameters k of the wire were 0.25, 0.4, and O.7 mm, 
respectively, and the distances x of the wire from the leading edge were 
15^ 30, 45, and 60 cm, respectively. When the wind speed V was low, the 
boundary layer was laminar all along the plate, but from a certain speed 
upward, the transition to turbulent flow was observed at that point where 
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the vire vas placed. Transition was detected hy a svidden change in the 
value of the total pressure G* in terms of the dynamic pressure q of 
the undisturbed stream, G* being indicated by a small pitot tube with 
a flattened mouth of 1 mm external width and O.3 mm width, which was 
placed in contact with the plate at a point 70 cm behind the leading 
edge. A sample record of measurements for k = 0.4 mm is shown in 
figure 29. 


When a flat plate is placed along a uniform stream of velocity V, 
the Blasius solution (ref. 40) of the laminar boundary layer equation 
gives 



for a point of distance x from the leading edge. Writing K for the 
critical value of kv^^v, the permissible height of projection k is 
given by 

0.576 1 = 

We determine V from the wind speed corresponding to the kink of the 
curve as shown in figure 29, plot 0.576(k/x) in a logarithmic scale 
against Vx/v, and draw a straight line of the elope -3/4 through the 
points (fig. 30) . We thus obtain K = I3, which is far greater than 
the value K = 7 estimated previously. 

32. Similar measurements were also performed on an airfoil section 
L.B. 24. The model was of 0.8 in span, of 1.2 m chord, fitted with end 
plates 1.3 ® X 0.6 m, and set at zero angle of attack in the same wind 
tunnel. Wires of various diameters (k = 0.25, 0.4, 0.7 mm) were attached 
parallel to the span, in contact with the svirface, at 10 percent of the 
chord from the leading edge (x = O.l). Iransition was detected by the 
sudden change in total pressure as indicated by a pitot tube with a 
flattened mouth of 2.7 mm external width and 0.9 mm depth, which was 
placed in contact with the surface at 50 percent of the chord from the 

leading edge. Results of measurements are shown in figiore 31, where 
V-f y 

R = and t is the chord length. 

The friction velocity may be generally expressed in the form 

V* = AVR“^/^ 
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where A is a function of s/t (s is the length measured along the 
surface from the fonfard stagnation point) . We can calculate A hy 
applying either the Pohlhausen approximate solution (ref. 4l) or the 
singplified method due to the author (refs. 15 and 16) . The permissible 
height of projection is then determined by 



Applying the Pohlhaiisen method to the theoretically calculated 
distribution of pressure, we get the values of function A, as shown 
in figure 32. Since A is 1.23 at the position of wire (x = O.l), 
the value 1.23(k/t) is plotted in a logarithmic scale against R in 
figure 33> R being the Reynolds number corresponding to the kink of. 
the curve as given in figure 31» It will be seen that although the 
measured points are on a strai^t line of the slope -3/^^ they give 
K = 15, which is somewhat hi^er than the value foimd for the flat plate. 

33* Now we apply the preceding result to the fragmentary data 
"hitherto known in order to check the adequacy of the estimate. First, 
we examine the results of wind tunnel experiment on a symmetrical 
laminar-flow airfoil L.B. 27>^ on which various projections are attached 
at 3 percent of the chord from the leading edge (x = O.O3) . The model 
was of 0.8 m span, of 1.2 m chord, and set at zero angle of attack. 

The profile drag was measvued by the method similar to that for L.B. 24 
(see paragraphs 12 to I3) . The res\ilts are shown in figure 34, from 
which we find that the rubber tape of O.O7 mm thickness gives no effect 
over the range of Reynolds nimabere R covered by the experiment, while 
the piano wire of 0.5 mm diameter gives a con5>letely turbvilent friction. 
The effect of the wire of 0.25 mm diameter begins to appear at 

R = 1.3 X 10^. Inserting the values k = 0.25 mm, t = 1.2 m and 

A = 1.95 in the formula A ^ = KR"^/^, we get K = 15.5> which is in 

good agreement with the result in paragraph 32. The value A = 1.95 
was read from figure 32, since the leading edge portion of L.B. 27 
almost coincides with that of L.B. 24. Profile drag coefficients of 
both airfoils are also the same over the range of Reynolds numbers 
examined. 

34. In connection with the determination of boundary layer transi- 
tion on airfoils in the MCA full-scale wind tunnel (ref. 42), an aux- 
iliary measurement has been reported, in which the effect was examined 


^L.B. 27 has a maximum thickness of 10 percent of the chord at 
60 percent of the chord from the leading edge. See figure 5 of paragraph 6, 




NACA TM 1351 


33 


of rubber tapes attached at 5 percent of the chord (x = O.O5) from the 
leading edge of the model. The airfoil section was NACA 0012, the chord 

was 72 inches, and the Reynolds nvimber R was 4.l8 x 10^. No effect 
was found when the tape was O.OO3 inch thick, some effect began to 
appear when it was O.OO6 inch thick, and the transition moved right to 
that position where the tape of O.OO9 inch thickness was attached. 
Assuming A = 1.6 and K = I5, we estimate from the preceding formula 
the value O.OO7 inch for the permissible thickness, which seems to agree 
well with the observation. 


If we further assume that the value A = 1.6 is also applicable 
to the case of the British fli^t experiment mentioned in paragraph 30, 
we find O.OOk inch for the critical height for transition with 

= 10 feet, R = 1,8 X 10*^, and K = 15. On the other hand, we esti- 

lev 

mate the permissible limit in turb\ilent boundary layer by * =4, 

which may be written in the form 


k 



by the relation ^ c^U]^^, where u^ is the local wind speed and 

is the coefficient of local skin friction. If we assime u-j^ = 1.2V 
and Cf = 0.003 (which is equal to the coefficient of mean skin friction 

for a flat plate at R = 1,8 x lo'^) , we obtain 0,0006 inch for the per- 
missible routes s thickness. Since the thickness of the camouflage 
psiiit is reported to be about 0,001 inch, it may be concluded that the 
paint Increases only the friction in turbulent boundary layer, without 
affecting, however, the transition to turbulent flow. This is in good 
agreement with the experimental resvilts. 
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APPENDIX 

TRANSITION CAUSED BI LAMINAR SEPARATION 


35« As Is well-known, the phenomenon of sudden decrease in drag of 

a sphere at a certain value of the Reynolds numher R = ~ (V is the 

speed of undisturbed stream and d is the diameter of sphere) is 
explained by supposing that the boundary layer separates while it is 
laminar when R is low, but it separates after transition to turbulence 
when R is hi^, thus restilting in diminishing the so-called dead water 
region. Probably the transformation from laminar separation to turbulent 
separation may proceed as follows: 


When the laminar boundary layer separates from the surface, the 
detached layer remains also laminar at first, but it is so unstable that 
it becomes turbulent at a short distance. This transition from laminar 
to turbulent flow is considered to occur when the local Reynolds number • 
based on the width of the detached layer and the velocity outside the 
layer exceeds a certain value, so that the transition moves upstream 
toward the separation point as R increases. When the transition 
approaches sufficiently near the separation point, it becomes possible 
for the detached layer to come back again to the downstreeim surface, 
because the turb\ilence produced will drive the flow forward. The layer 
reattaches to the surface as a turbulent layer, and accordingly the drag 
coefficient begins to decrease. The distance between the separation and 
the first turbulent boundary layer decreases as R increases, and finally 
the fully developed turbulent boundary layer commences just downstream 
of the separation point. The drag coefficient then ceases to decrease. 


36. Now, in order that the separated layer reattach to the surface, 

U-10 

it seems necessary for the local Reynolds number Rq = at separation 

to exceed a certain critical value, where uq^ is the velocity outside 

the boundary layer, and 9 is the momentum thickness of the boundary 
layer. This may be explained as follows: According to the laminar 

boundary layer theory, the separation occurs when the quantity 


^ ds <is ’ 9 


exceeds a certain value, suggesting that the pressure rise (dp/ds)0 . 
becomes too large in proportion to the shearing stress at the sur-? 
face pv(ui/0). Assuming analogically that the separated layer leaves 
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the surface when the pressure rise becomes too large in proportion to 

the momentum pu^ , we then find that Rg at the separation point must 

exceed a certain critical value in order that the separated layer 
reattach to the surface. 

37 . In order to determine the critical value of Rq, we consider 

in detail the condition where the coefficient of sphere drag begins to ' 
decrease. This condition corresponds to the point B of the curve of 
figxore 35^ which represents an idealized variation of the drag coeffi- 
cient Cj)j^ drag j- (pV2/2j ^nd^A^ or the pressure difference coeffi- 
cient Ap/q|= difference of pressures at the forward stagnation point 

and the point corresponding to the central angle 157 . 5 ° — (pV^/2^ with 

the Reynolds number R. Within the range AB, the pressure distribution 
aro\md the sphere is approximately independent of R; the typical example 
may be found from the experiments due to Fage (ref. 43 ). Fortunately, 
the boxmdary layer calculation has also been performed for that distri- 
bution of pressiare by Tomotika and Imai (ref. 44), so that the local 
Reynolds number Rq is given by 


Re = o.4o\|r 


at the separation point. Although the calculation has originally been 

made for a particular Reynolds number, R = I .57 x 105, the preceding 
relation may be applied for any value of R in the range of AB. Putting 
the value of R at B, and writing 


% 


crit 



we have as the critical value of Rg above which the separated 

laminar layer reattaches to the surface. Conventionally the Reynolds 
number Rc corresponding to Cjj = 0.3 (or Ap/q = 1.22) has been used, 

instead of R^, for representing the degree of stream turbulence, but it 
is not so difficult to estimate the value of Rg from the measured cxirve 

Ap/q) against R. For example, we have from the experiments 
of towing spheres in the free atmosphere (ref. 35 ) 

, Rb = 3.6 X 1q 5 Rp = 3.85 X 105 

®crit 


= 240 
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Since these values refer to the case of very low turbulence, it will be 

seen that Rg . = 2l(0 represents the hipest possible value. It is 

crit 


also to be noted that 


He 


0.9^ in this example and that almost the 


same value has been obtained by the author's experiments on spheres of 
various diameters (ref. 29) . 


38. If the stream turbulence is not low so that Rq is less than 
3.85 X 10^, then will be less than 2^. Assuming the ratio Rg/Rc 

to be constant, we can estimate the corresponding value by 


Rfl = 240| 
^crit 


R. 




10 ' 


On the other hand, we can also estimate the value of H^crit <ilr®ctly 

from the boundary layer measurements. The results of the analysis for a 
sphere as well as circular and elliptic cylinders are simnnarlzed in the 
following table, where is the critical value Rq estimated 

from Rp by the preceding formula, R^ is the value of Rg observed 

sep ^ 

at the separation point when the boundary layer really separates while 
it is laminar, and Rq is the value of Rg observed at the calcu- 

lated laminar separation point when the boundary layer separates after 

transition. The fact that Ra .. lies between Ra and Ra. 

''crix v^sep trans 

seems to suggest the adequacy of the preceding consideration. 


Body 

He 

^crit 

Re 

^sep 

^ trans 

Reference 

Sphere 

2.5 X lo5 

190 

160 

220 

43 

Circular cylinder 

1.5 X 1 q 5 

150 

i4o 

225 

1^5 

Elliptical cylinder 

2.7 X io5 

200 

160 

koo 

^6, 47 


39* now proceed to apply our result to interpreting the effect of 
Reynolds number on maxlmtmi lift of airfoils. For the angle of attack 
near the stall, the flow separates shortly downstream of the leading edge 
while the boundary layer is laminar. If the flow fails to reattach to 
the surface as a turbulent layer, the maximum lift coefficient of 
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the airfoil will be almost independent of the Reynolds number R, At 
the separation point, similar to the case of a sphere, the relation of 
the form 


Rq = klfR 


holds, where R = ^ is the Reynolds number referred to chord length t,® 

and k is a constant depending on the shape of airfoil and the value of 
^ that Rq is less than the critical value Rg 

crit 

then ^LniftY independent of R. Assuming the same value of 

Rq . . for the sphere as for the airfoil, we obtain 
crit ' 


R = 


0.16 

k2 


% 


for the Reynolds nmber above which begins to increase with the 

Reynolds number. Therefore, the ratio of the Reynolds number corresponding 
to a certain value of of an airfoil and the critical Reynolds 

number of sphere in the same stream, Rg or Rq , becomes independent of 

the stream tvirbulence. Denoting the values for a reference tunnel with 
asterisk, we have 


R Rc 


which in iurn means “thai ihe ratio of Reynolds numhers corresponding to 
a certain value of is equal to the ratio of critical Reynolds 

nvunbers of a sphere. This is useful for comparing the values of 

obtained in two different wind tunnels. Considering the reference condi- 
tion to be the free flight in the atmosphere, we find 

R* = R X ^ 

Rc 


®It is to be noted that R is referred to t, while Rg and Rq 
are referred to d. 
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as the free, flight Reynolds number vhich will give the same value of 
*^^aax 'that observed in a wind tunnel. .This is just what is called 

the effective Reynolds number. Strictly speaking, such an argument as 
mentioned before should apply only to the Reynolds numbers near the 
critical value, but there are many experimental evidences showing the 
usefulness of the concept of effective Reynolds number for most practi- 
cal purposes, as far as the commonly used airfoils and range of Reynolds 
numbers of both wind tunnel and free flight are concerned. 

40. Finally, we consider a more quantitative example to show the. 
adeqxiacy of the preceding argument. In figure 36, various 

MCA symmetrical airfoils are plotted against the effective Reynolds 
number R*, the experimental data being taken from the results of the 
MCA variable-density wind tunnel (ref. 48). As already mentioned, iq> 
to a certain value of R*, almost independent of R*» This 

corresponds to the condition in which the laminar separation just behind 
the leading edge fails to reattach to the surface, resulting in a con- 
siderable dead water region above the airfoil surface. The value of 
Chnifly approximately 0,9, irrespective of the thickness; it is almost 

equal to the value for a flat plate of vanishing thickness. Theory of 
discontinous flow, when applied to the flat plate, seems to give a lift 
coefficient close to O.9 (ref. 4-9). We idealize, therefore, the experi- 
mental curve as shown by dotted lines in figure 38. Then, the point 
where the dotted line meets the line C^^^^ = 0,9 will be considered to 

correspond to Rfl = 240. In order to determine this point, we calcu- 
crit 

late the value of Rq at the laminar separation point for a lift 
coefficient Cl = O. 9, and the value of R which gives Rg = 240. We 
first calculate the pressure distribution by the formula of paragraph 3 
for the two-dimensional potential flow around the airfoil section. 
Althou^ the formula may be applied to any arbitrary airfoil section, 
we have determined the pressure distribution only for the airfoil 


T = 0.287erfx(l - x)(5 - 4x) 


in order to sin5)lify the calcvilation, because no great exactitude is 
required in the present problem, x is the coordinate along the chord, 

<x = 0 »ni^ X = 1 corresponding to the leading and trailing edges, 
respectively, T is the half-thickness, and e is the maximum thick- 
ness in terms of chord length. The airfoil represented by the preceding 
expression coincides with sufficient accuracy with the true MCA symmetri- 
cal airfoil except near the trailing edge. The value of Re/|^ at the 
laminar separation point was then determined for the calctilated pressure 
distribution by applying the approximate method due to the author (refs. 15 
and 16) . 
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The critical values 'thus calculated are shovn by a solid 

line in figure 37^ while the corresponding values taken from figure 36 
are shown by white circles. The agreement is fairly good, and especially 
satisfactory when the thickness of the airfoil is small. In general, the 
thickness of the boundary layer near the trailing edge increases as 

increases. If is further increased, however, a laminar separation 

suddenly occurs near the leading edge when the thickness is small, while 
the trailing-edge turbulent separation moves a considerable extent for- 
ward before the leading-edge laminar separation occurs when the thickness 
is moderate. Therefore, the assumption of the analysis is more satis- 
factorily realized in the case of small thickness, thus bringing the 
calciiLated and observed values in close agreement. 

In conclusion, the author wishes to acknowledge his indebtedness 
for the assistance given by Messrs. C. Noda, S. Mituisi, I. Shinra, 

S. Asaka, R. Hama, and K. Takeda. 


Translation by Itiro Tani 
University of Tokyo 
Tokyo, Japan 
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TABLE 1.- ORDINATES FOR SYMMETRICAL SECTION. VALUES OF T/e 


X 

I 


K 

L 

M 

N 

0 

0 

0 

0 

0 

0 

0 

,003 

.0465 

.0563 

.0574 

.0584 

.0606 

.0627 

.006 

.0661 

.0793 

.0808 

.0823 

.0854 

.0886 

.0125 

.0962 

.1135 

.1158 

.1180 

.1226 

.1275 

.025 

.1374 

.1-589 

.1620 

.1652 

.1720 

.179^ 

.05 

.1963 

.2208 

.2252 

.2298 

.2396 

.2508 

.075 

.2414 

.2663 

.2716 

.2771 

.2892 

.3031 

.10 

.2789 

.3029 

.3089 

.3151 

.3288 

.3447 

.15 

.339*^ 

.3599 

.3667 

.3738 

.3893 

.4072 

.20 

.3865 

.4026 

.4097 

.4170 

.4328 

.4503 

.25 

.4236 

.4353 

.4422 

.4492 

.4638 

.4788 

.30 

.4523 

.4601 

.4663 

.4724 

.4844 

.4949 

.35 

.4737 

.4783 

.4833 

.4881 

.4962 

.5000 

.40 

.4885 

.4906 

.4942 

.4971 

.5000 

.4953 

.45 

.4972 

.4977 

.4994 

.5000 

.4948 

.4818 

.50 

.5000 

.5000 

.4994 

.4931 

.4797 

.4604 

.55 

.4970 

.4944 

.4873 

.4781 

.4558 

.4320 

.60 

.4871 

.4778 

.4656 

.4522 

.4241, 

.3975 

.65 

.4691 

.4509 

.4342 

.4175 

.3856 

.3578 

•.70 

.4418 

.4139 

'.3939 

.3751 

.3413 

.3138 

.75 

.4038 

.3675 

.31^55 

.3258 

.2921 

.2664 

.80 

.3538 

.3121 

.2899 

.2706 

.2393 

.2165 

.85 

.2908 

.2481 

.2277 

.2105 

.1836 

.1650 

.90 

.2133 

.1762 

.1598 

.1464 

.1262 

.1128 

.95 

.1201 

.0966 

.0870 

.0793 

.0680 

.0609 

1.00 

.0100 

.0100 

.0100 

.0100 

.0100 

.0100 
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TABLE 2.-' AUXILIARY FUNCTION, Bg ASSOCIATED WITH THE 
. PRESSURE DISTRIBUTION OF SYMMETRICAL SECTIONS 


X 

I 

j 

K 

L 

M 

N 

0.0125 

0.59 

1.10 

1.13 

1.16 

1.20 

1.24 

.025 

.66 

1.06 

1.09 

1.13 

1.20 

1.26 

.05 

.75 

1.04 

1.07 

1.11 

1.20 

1.29 

.10 


1.03 

1.07 

1.11 

1.21 

1.33 

.20 

.96 

1.03 

1.08 

1.12 

1.22 

1.36 

.30 

1.02 

1.04 

1.09 

1.14 

1.23 

1.33 

.40 

1.05 

1.05 

.1.10 

1.15 

1.23 

1.22 

.50 

1.08 

1.17 

1.19 

1.19 

1.14 

1.04 

.60 

1.12 

1.18 

1.11 

1.05 

.92 

.80 

.70 

1.09 

.98 

.88 

.79 

.63 

.51 

.80 

.89 

.65 

.52 

.42 ♦ 

.28 

.19 

.90 

.33 

.06 

-.03 

-.10 

-.18 

-.22 

.95 

-.32 

-.47 

-.50 

-.52 

-.53 

-.52 

.975 

-.95 

-.93 

-.91 

-.89 

-.84 

-.78 




NACA TM 1351 


TABLE 3.- AUXILIARY FUNCTION Be ASSOCIATED WITH THE 
PRESSURE DISTRIBUTION OF SYMMETRICAL SECTIONS 


X 

■ I 

J 

K 

L 

M 

N 

0.0125 

-L .04 

- 4.40 

- 4.88 

-4.90 

-5.20 

-5.30 

.025 

-2.86 

-3.08 

-3.20 

-3.22 

-3.37 

-3.55 

.05 

- 2.04 

-2.08 

-2.13 

-2.17 

-2.28 

-2.39 

.10 

-1.39 

-1.33 

-1.35 

-1.37 

- 1.42 

- 1.48 

.20 

-.83 

-.74 

-.74 

-.74 

-.74 

-.70 

.30 

-.50 

-.43 

-.41 

-.39 

-.31 

-.21 

.1^0 

-.24 

-.19 

-.17 

-.12 

.00 

.18 

.50 

.00 

.00 

.11 

.22 

.39 

. .50 

.60 

.27 

.44 

.53 

.60 

.70 

.75 

.70 

.65 

.83 

.89 

.92 


.91 

.80 

1 . 11 , 

1.19 

1.18 

1.15 

1.09 

1.02 

.90 

1.70 

1.52 

l. 4 l 

1.31 

1.16 

1.04 

.95 

2.03 

1.66 

1.50 

1.36 

1.16 

1.03 

.975 

2.20 

1.73 

1.54 

1.38 

1.15 

1.02 
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TABLE 4.- ORDINATES FOR MEAN CAMBER LINES. VALUES OF M/f 


X 

1^0 


®3 

1^5 


Ej^ 

^3 

0 

.003 

.006 

.0125 

.025 

.050 

.075 

.10 

.15 

.20 

.25 

.30 

.35 

.40 

.45 

.50 

.55 

.60 

.65 

.70 

.75 

.80 

.85 

.90 

.95 

1.00 

X for M/f = 1 

0 

.0295 

.0529 

.0969 

.1687 

.28a 

.3843 

.4690 

.6098 

.7219 

.8113 

.8813 

.9341 

.9710 

.9928 

1.0000 

.9928 

.9710 

.9341 

.8813 

.8113 

.7219 

.6098 

.4690 

.28a 

0 

.500 

0 

.0282 

.0509 

.0914 

.1636 

.2797 

.3772 

.4620 

.6045 

.7191 

.8112 

.8838 

.9384 

.9758 

.9962 

.9988 

.9813 

.9414 

.8794 

.7878 

.6898 

.5673 

.4305 

. 2^5 

.1366 

0 

.482 

0 

.0292 

.0526 

.0969 

.1695 

.2897 

.3907 

.47a 

.6252 

.7424 

.8353 

.9067 

.9579 

.9892 

1.0000 

.9881 

.9472 

.8735 

.7704 

.a 59 

.5108 

.3766 

.2535 

.1493 

.0651 

0 

.450 

0 

.0300 

.oao 

.0995 

.1740 

.2973 

.4006 

.4902 

.6396 

.7579 

.8507 

.9207 

.9688 

.9953 

.9988 

.9761 

.9178 

.8192 

.6903 

.a 74 

.4078 

. 2 a 7 

.ia 6 

.1071 

.0472 

0 

.433 

0 

.0371 

.0666 

.1220 

.2119 

.3585 

.4787 

.5808 

.7443 

.8638 

.9451 

.9898 

.9963 

.9589 

.8625 

.6300 

.3885 

. 2 a 6 

.1811 

.1207 

.0772 

.0459 

.0242 

.0101 

.0024 

0 

.333 

0 

.0313 

.05a 

.1039 

.1814 

.3093 

. 4 i 6 i 

.5084 

.6615 

.7812 

.8733 

.9387 

.9826 

.9998 

.9872 

.9387 

.8561 

.7473 

.6259 

.4918 

.3682 

.2590 

.1686 

.0973 

.0425 

0 

.406 

0 

.0327 

.0589 

.1081 

.1887 

.3212 

.4315 

. 5266 
.6831 
. 804 l 
.8952 
.9586 

.9939 

.9973 

.9679 

.9085 

.8238 

.7204 

.6056 

.4871 

.3721 

.2669 

.1761 

.1020 

.0442 

0 

.381 
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TABLE 5.- AUXILIARY FUNCTION Ac ASSOCIATED WITH THE 
PRESSURE DISTRIBUTION OF MEAN CAMBER LINES 


X 

Do 

Dl 

D3 

D5 

D.0 

% 

^3 

0.0125 


-0.53 

-2.50 

-3.17 

-5.49 

-3.51 

-3.62 

.025 

*^.53 

.88 

-.47 

-.91 

-2.17 

-1.07 

-1:08 

.05 

^^.53 

1.89 

.99 

.71 

.22 

.67 

.74 

.10 

1^.53 

2.61 

2.03 

1.88 

1.94 

1.93 

2.05 

.20 

^.53 

3.15 

2.80 

2.75 

3.20 

2.85 

3.01 

.30 

1^.53 

3.40 

3.17 

3.15 

3.79 

3.29 

3.46 

.ItO 

4.53 

3.56 

3.40 

3.41 

4.18 

3.57 

3.61 

.50 

4.53 

3.68 

3.57 

3.61 

4.46 

3.50 

3.10 

.60 

4.53 

3.61 

3.21 

2.95 

-1.03 

2.35 

2.14 

.70 

4.53 

3.19 

2.07 

1.31 

-.83 

1.00 

1.04 

.80 

4.53 

2.43 

.75 

.05 

-.63 

.03 

.17 

.90 

4.53 

1.34 

-.05 

-.26 

-.42 

-.26 

-.22 

.90 

4.53 

,68 

-.15 

-.20 

-.29 

-.21 

-.21 

.975 

4.53 

.33 

-.08 

-.14 

• -.20 

-.15 

-.15 






TABLE 6.- AUXILIARY FUNCTION Ag ASSOCIATED WITH THE 
PRESSURE DISTRIBUTION OF MEAN CAMBER LINES 


X 

Do 

Dl 

l >3 



E]^ 

^3 

0.0125 

-6.31 

-6.13 

-6.35 

-6.52 

-7.92 

- 6,19 

-7.07 

.025 

-5.28 

-5.18 

-5.37 

-5.50 

-6.61 

-5.73 

-5.94 

.05 

- 4.25 

- 4.21 

-4.36 

- 4.46 

-5.25 

-4.63 

- 4.79 

.10 

-3.17 

-3.19 

-3.29 

-3.36 

-3.78 

- 3.46 

-3.55 

.20 

-1.99 

-2.05 

-2.09 

-2.10 

-2.00 

-2.11 

-2.11 

.30 

-1.21 

-1.27 

-1.22 

-1.18 

-.52 

1.09 

-.99 

.Uo 

-.58 

-.58 

-.43 

-.31 

1.26 

-.07 

.26 

.50 

0 

.14 

.50 

.76 

00 

1.33 

1.46 

.60 

.58 

1.02 

1.79 

2.32 

2.00 

2.38 

2.21 

.70 

1.21 

1.90 

2.63 

2.87 

1.02 

2.30 

2.36 

.80 

1.99 

2.61 

2.60 

2.24 

.52 

2.00 

1.97 

.90 

3.17 

2.97 

1.88 

1.35 

.21 

1.25 

1.31 

.95 

4.25 

2.91 

1.63 

1.05 

.10 

.91 

1.01 

.975 

5.28 

2.75 

1.59 

.94 

.05 

.82 

.88 




















TABLE 8.- AIRFOIL WITH UNIFORM DISTRIBUTION OF PRESSURE 


NACA TM 1351. 


cnt-\c\o ^ 

ONco t^Crvcvjoifxou^ooo 

0 \ On ON 0 \ On CaOO CO t^VD uS 
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-4- 00 

o 00 on o CO o cvj coco ^ iH o 

• • On 00 t — t**- 

fHrH •••••••••• 


5>jSI «HpH^c^OMrNonHoo 
CO 00 iH cn on on on m on 

o ON On CO ooScocococococo 


on on ON iPv utn^ -it co m 
ONlfNoncVirHHHHHIHHH 
0\0\<y\<T\0\0\0\0\0\0\0s0\ 


On ONrHVOO®^ CVJO o 
COVO^OJCVJf-IHOOOOO 

•••••••••••A 


iH irNoo-=t itvonvo cvjno \r\^ a 
^CO CVJ0-^0'^-:tCVJrH00 
VO-:tcnCVjrHHOOOOOO 

• •• ••••••••• 

Hi— I fHiHiHHrHHrHHrHrH 


ON>-CO O^-lTNoncVjHOOO 
onOJH HOOO 00000 

• ••••«•••••• 

HHHHHHHHHHHH 


COt^lTNOlTNOOOOOO 
C?N ON ON CTnOO CO C^VD OTS^ CVJ 


VO CVJ ITN O ^ 

OHCVjLTNC^OirNol^OOO 
OOOOOHHcCiCVJ fO-^ uTN 
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Figure 9 
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Figure 18 
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Figure 20 







Figure 22 
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